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Abstract 

We investigate the thermodynamic curvature resulting from a Riemannian ge- 
ometry approach to thermodynamics for the Pauli paramagnetic gas which is a 
system of identical fermions each with spin i. We observe that the absolute value 
of thermodynamic curvature can be interpreted as a measure of the stability of the 
considered system only in the classical and semiclassical regime. But in quantum 
regime some exceptions are observed. 

1 Introduction 

Thermodynamic fluctuation theory whose basic goal is to express the time independent 
probability distribution for the state of a fluctuating system in terms of thermodynamic 
quantities, is usually attributed to Einstein who applied it to the problem of blackbody 
radiation [1]. The full formalism for classical thermodynamic fluctuation theory was 
worked out by Green and Callen [2] in 1951 and elaborated upon by Callen [3]. 

However, despite of a wide range of applicability, the classical fluctuation theory fails 
near critical points and at volumes of the order of the correlation volume and less. 

In 1979 Ruppeiner [4] introduced a Riemannian metric structure representing ther- 
modynamic fluctuation theory, and related to the second derivatives of the entropy. His 
theory offered a good meaning for the distance between thermodynamic states. He showed 
that the breakdown of the classical theory occurs because it does not take into account 
local correlations [5]. This deficiency of the classical theory is precluded in the covariant 
fluctuation theory of Ruppeiner by using a hierarchy of concentric subsystems, each of 
which samples only the thermodynamic state of the subsystem immediately larger than 
it [6,7]. One of the most significant topics of this theory is the introduction of the Rieman- 
nian thermodynamic curvature as a qualitatively new tool for the study of fluctuation 
phenomena. It is this geometry which is the basis for this paper. Here we investigate the 
case of a Pauli paramagnetic gas that is a gas of identical spin | fermions in the presence 
of an external magnetic field. 
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The outline of this paper is as follows. First the Riemannian geometry of thermody- 
namic fluctuation theory is summarized. Second, the Riemannian scalar curvature of the 
Pauli paramagnetic gas is evaluated. Finally the curvature of the classical ideal param- 
agnetic gas is calculated and is compared to the the curvature of Pauli paramagnetic gas 
in the classical limit. 



2 Geometrical view point of thermodynamics 

In this section we review the Riemannian geometry of thermodynamics, discuss its con- 
nection to the covariant thermodynamic fluctuation theory, and summarize the resulting 
interpretation of thermodynamic curvature. Riemannian structure of the thermodynamic 
state space is defined by the second derivatives of a thermodynamic potential density as a 
metric tensor [7]. If we choose extensive densities as coordinates, we can use either energy 
or entropy density as the potential and these two descriptions are thermodynamically 
equivalent; but the metric tensor is different in these two representations. Here we work 
in entropy representation, because then the meaning of the distance, measured in units 
of average fluctuations, is very transparent. When some extensives are substituted by 
intensives, the potential is a Massieu function [3]. Now we consider an open subsystem 
Ay, with fixed volume V, of a thermodynamic fluid system Ay with a very large volume 
Vq. The system Ay consists of r fluid components and is in equilibrium. We denote by 
the n-tuple ao = (a°, aj, a 2 ,, • • • , ciq) the internal energy per volume and the number of 
particles per volume of the r components of A„ [7]. These parameters are the standard 
densities in the entropy representation; they constitute the thermodynamic state of Ay . 
The subsystem Ay has the corresponding thermodynamic state a. The Gaussian approx- 
imation of the classical thermodynamic fluctuation theory asserts that the probability of 
finding the thermodynamic state of Ay between a and a + da is[7]: 

P v (a\a )da°da 1 ■ ■ ■ da r = (^—^j exp[-—g^ u (a )Aa' l Aa' / } 

x ^g(a )da°da 1 ■ ■ ■ da r , (2.1) 

where 

Aa^ = a^-a^ 



1 d 2 s 
9>1U ~ Y^da^da" 



(2.2) 



a=ao 



where s is the entropy per volume in the thermodynamic limit, Kb is Boltzman's constant, 
and g(a ) = det[^(a )]. 

The quadratic form in Eq.(2.1), 

(A/) 2 = ^(ao)A^A^ (2.3) 

constitutes a positive definite Riemannian metric on the thermodynamic state space. 
The positive definiteness results since the entropy is a maximum in equilibrium a = 
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a . Eqs.(2.1) and (2.3) denote the Physical interpretation for the distance between two 
thermodynamic state. The less the probability of a fluctuation between the states, the 
further apart they are. The quantity 




'^da 1 ■■■da r 



in Eq.(2.1), is the invariant Riemannian thermodynamic state space volume element. The 
form of Eq.(2.2) holds only in standard densities. To express the metric tensor in a general 
set of thermodynamic coordinates x — x(a), one can use the following transformation rule 



Having the metric we can calculate the Riemannian curvature tensor. For our metric, the 
scalar curvature R ,has units of real space volume, regardless of the dimension of the state 
space[ll]. It is a measure of effective interaction between the components of the system, 
proportional to the correlation volume, and diverges near the critical point of the pure 
interacting fluid. Covariant thermodynamic fluctuation theory indicates that curvature is 
a measure of the smallest volume where classical thermodynamic fluctuation theory could 
work. This theory was proposed as the correct way to extent the classical thermodynamic 
fluctuation theory beyond the Gaussian approximation [7]. 

An alternative interpretation of the thermodynamic curvature was offered by Janyszek 
and Mrugala[8]. They suggested that the thermodynamic curvature is a measure of the 
stability of the considered system. The system is less stable if the curvature increases and 
vice versa. Also these authors calculated the curvature of ideal Fermi and Bose gases [9]. 
They show that these systems have the curvature with opposite signs. 

In this paper we interpret the absolute value of the curvature as a measure of stability 
in order to come to an agreement with the curvature of the boson ideal gases that diverge to 
negative infinity (in the sign convention used here) where the Bose-Einstein condensation 
occurs [9]. 



We now turn our attention to studying the equilibrium state of a gas of noninteracting 
fermions in the presence of an external magnetic field H. 

The extensive parameter which describes the magnetic properties of a system is M, 
that is the component of the total magnetic moment parallel to the external field. The 
entropic intensive parameters are defined as [3], 




<9a M da' 



(2.4) 
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(3.1) 



We use the thermodynamic potential <f> which is defined as, 



where u,p,m and P are energy per volume, density, magnetization and pressure respec- 
tively. The energy of a particle, in presence of an external magnetic field H, is given 
by 

c P -J.H (3.3) 



2m, 

where J is the intrinsic magnetic moment of the particle and mo is its mass. For the case 
of the Pauli paramagnetic gas the spin of each particle is |; the vector J must then be 
either parallel to the vector H or anti parallel. From the grand canonical distribution 
(using Fermi-Dirac statistics) one can obtain the following equations [13]: 

inQ = < 3 - 4 > 



2 



where 



./;; = /nOf); ( (3.6) 

+ r l r A 4 1 I rr\ 

V = ^xp[ T — ] = exp[— T —) (3.7) 

and A = r is the mean thermal wavelength of the particle, h is the Planck 

{2nm K B T)? 

constant and 

we have used the standard symbol for the fugacity rj = exp(-^A^). From Eqs.(3.2) and 
(3.4) the thermodynamic potential is obtained 

<f>(x,y,z) = /x-S[/ f (e-^ J *) + fi(e-y +Jz )}. (3.9) 

Where / = and x = F l ,y = F 2 , z = F 3 . We have set K B = 1. 

Now it is straightforward to obtain the metric elements in F coordinates [7]: 

d 2 4> 

9»» = avuaw ( 3 - 10 ) 



dF»dF v 
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drj 



According to Eqs.(3.9) and (3.10) and noting that, = ^fn-i(v)i the components of 

the metric tensor are as follows 

9n = ^-Ix~i(a + b) g 12 = ^Ix~^(c + d) g 13 = --IJx~*(c- d) 

922 = Ix-l(e + f) g 23 = -IJx-l(e-f) g 33 = IJ 2 x'ke + /) (3.11) 

where a = ft, b = /J, c = ft, d = fi, e = ft, f = fl; they are functions of y and z. 

2 2 2 2 2 2 

Their derivatives with respect to y and z are as follows: 

|2 = -c ^ = Jc ^ = -d f = -3d 

ay oz ay oz 

fy = -C f,=J* fy = ~f fz=~Jf (3-12) 

oy oz oy oz 



where h = f\, k = these quantities are used for obtaining the scalar curvature. 

2 2 

The Riemann and the Ricci tensors and the scalar curvature are respectively 

R = g^Ric^ (3.13) 

where Ts are the Christoffel symbols[14]. The scalar curvature may be worked out with 
Eqs.(3.11), (3.12), (3.13): 

R = 2(5 e / a + 5e/ 6 -3eV-3^ e) ^ |55/V + 55/2fe2 

- 28 fee 2 - 25f 2 ach - 25f 2 bch - 28fd 2 e 2 + 12fd 2 ch 

- 25e 2 dka - 25e 2 bdk + \2c 2 dke + Ihcdhka + Ihcdbhk] (3. 14) 

As it is seen from Eq.(3.14), R is a symmetric function of z; this means that scalar 
curvature is independent of the orientation of external magnetic field, R(—H) = R(H). 

In the classical limit in the absence of the external magnetic field, we have r; ± — > 77 — > 
and fniv) ~~ * V'i so R * s obtained as follows: 

1 A 3 

R=-— (3.15) 
4 77 



On the other hand in this limit, Eq.(3.5) results 

2 



r v ,'/ (3-16) 

From Eqs.(3.15) and (3.16) the classical limit of R is given by 

R=± (3.17) 

This surprising simple result, shows that in the classical limit the scalar curvature is in 
the order of the volume occupied by a single particle. It is in complete agreement with 
the scalar curvature obtained by Ruppeiner for multicomponent ideal gas [11]. It means 
that in the classical limit the scalar curvature of the Pauli paramagnetic gas behaves like 
that of a two-component ideal gas. 

In Fig.l we present the dependence of R on rj for a fixed value of H, for an isotherm 
in units of A 3 . In the classical region where 77 < 1, R diverges near r] = 0. It relates to the 
fact that in this limit p goes to zero (as Eq.(3.17) demonstrates), so there are not enough 
particles for a continuous thermodynamic description. In the quantum mechanical region, 
where r\ 3> 1, R tends to a constant value. 

Fig. 2 shows the dependence of R on H for a fixed value of 77, for an isotherm in 
units of A 3 . R is a monotonically decreasing function of H. Physically, as the external 
magnetic field increases, the magnetization fluctuation decreases. So the system becomes 
more stable. Here we can interpret R as a measure of the stability of the thermodynamic 
system. The less the magnitude of R, the more stable the system becomes. 
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Fig.2 



In Fig. 3 there is something new that we can not explain it in terms of fluctuations. 
This figure shows the dependence of R on H for r\ = 10. As it is seen in this quantum 
regime, R has maximum around z — 2. Fig. 4 shows that this maximum of R occurs at 
higher values of z as we increase the value of rj. This behavior is in contrast with our 
intuition; because the magnetization fluctuation, < [^^ 2 > is a monotonically decreasing 
function of H for all values of r\. So it seems that the stability interpretation of thermo- 
dynamic curvature fails in strong quantum regime. 
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Fig.3 



Fig.4 



For the last point we allude to a relationship between R and the correlation volume. We 
note that the correlation function of Fermi gas in classical regime (pA 3 <C 1 or r\ <C 1) is 
given by the formula[12] 

(3.18) 



1 - 2 ~*J 2 

v{r) = — e ^ . 
v ; 2 



So one can see that the correlation volume in this regime is 

A 3 



V„ 



(27T 



(3.19) 



It means that in the classical regime Rr\ is proportional to the correlation volume. The 
possible relationship between the curvature and the correlation volume in the quantum 
regime has not been explored. 
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4 Classical ideal paramagnetic gas 



In this section we calculate the scalar curvature of a classical paramagnetic gas and 
compare it to the scalar curvature of the Pauli paramagnetic gas in the limit of low 
fugacity and low magnetic fields. 

Consider a gas of identical, mutually noninteracting and freely orientable dipoles, each 
having a magnetic moment J. In the presence of an external magnetic field H, the dipoles 
experience a torque tending to align them in the direction of the field. The energy of a 
particle is given by 

2 

r -JH cos (4.1) 



2m 

Here we have neglected the effect of the induced magnetic field. Mijatovic et al used the 
energy form of the metric to evaluate the geometry in the paramagnetic ideal gas[10]. 
They also used the particle number as the fixed scale. Here we use the energy form of 
the metric, we also use the volume as fixed scale. The logarithm of the grand canonical 
partition function (using Maxwell-Boltzman statistics) is obtained as follows: 

PV V smh( Jz) , AnS 

In Q c = = Ann— — - (4.2) 

Vc K B T 'A 3 Jz V ; 

and the thermodynamic potential is again obtained from Eqs.(4.2) and (3.2) 

3 „sinh( Jz) , , „ s 

C = A-Klx~ie- y j — '- (4.3) 

Where x,y and z are those that are defined in section 3 and Kb = 1. From Eq.(3.10) the 
metric elements are obtained 

1K T -i _„sinh(Jz) 
gn = 15nlx 2 e y — - 



Jz 

_s sinh(Jz) 

g 12 = 6nlx 2 e v ^ — l 

Jz 

,,cosh(Jz) sinh(J^). 
9is = -67r/x-3e-»(— ^ - — j^) 

. T _s sinh(Jz) 

Jz 

3 _„ .cosh(Jz) sinh(j2;). 

# 23 = -AtxIx 5 e y( 1 — L J^- 1 ) 
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3 „, J smh(Jz) cosh(Jz 
A-Klx~e- y ( i — - - 2 \ — 



+ (4.4) 
Using Eqs.(3.13) and (4.4), one can calculate the scalar curvature 

Rc = J_^_Il_ (4.5) 
8n r] sinh Jz 

Eqs. (4.3) and (4.5) show clearly that R c and 4> c satisfy the following equation. 

R c = ^ (4.6) 

<Pc 
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where k — \ and Kb = 1. This interesting result is nothing except the geometrical 
equation with n — |[7]. 

On the other hand , the equation of state of the classical ideal gas (PV = NK B T) 
and Eqs.(4.6) and (3.2) results 

R < = Y P < 47 > 

Eq.(4.7) shows that the curvature of the classical ideal paramagnetic gas is in order of the 
volume occupied by a single particle. 

We can see the magnetic field dependence of R c through Eq.(4.5). It is a monotonically 
decreasing function of z and has a maximum at z — 0. Here we can interpret the curvature 
as a measure of stability, since the magnetization fluctuation, < ^^a > , is dominant near 
z = 0, and decreases monotonically as z increases. 

Now let us look at the Eq.(3.14) for the curvature of the Pauli gas in the limit of low 
fugacity and low magnetic field, where — > ^ then one can obtain 

Eq.(4.8) is similar to Eq.(4.5), because the behavior of coshz is similar to that of 
but they are somewhat different. The source of this difference is related to the fact that 
in the case of the classical ideal paramagnetic gas, each dipole is freely orientable whereas 
each particle in the Pauli paramagnetic gas can choose only two directions (even in the 
limit of low fugacity) . 

In fact, had we constrained the dipoles to choose only two directions (parallel or anti 
parallel to H ) i.e. set cos6> = ±1 in Eq.(4.1) and used the classical grand partition 
function we would have obtained 

V 

In Q' c = 2r?— cosh Jz. (4.9) 



This is just Eq.(3.4) in the limit of low fugacity and low magnetic fields (where 



^ihy diiu. iuw iiiag,±icijiL. nciu.o ^vvncic j 

rf 1 = r]e ±zj ). The thermodynamic potential is 

<p' c = 2Ix^e~ y cosh Jz. (4.10) 

Again this is the classical limit of Eq.(3.9). Using Eqs.(3.10) and (3.13) one can obtain 
the scalar curvature which is just Eq.(4.8) 



5 Conclusion 

We have evaluated the thermodynamic curvature for the Pauli paramagnetic gas that is 
a system consisting of identical spin | fermions. 

In the classical limit (i.e 77 <§; 1) in the absence of the external magnetic field this 
curvature reduces to that of a two-component ideal gas. In this regime we can find a 
simple relationship between the curvature and the correlation volume. 
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In the limit of low fugacity and for a finite value of external magnetic field the curvature 
of the Pauli gas coincides with that of the classical ideal paramagnetic gas which is a 
monotonically decreasing function of magnetic field; here we can interpret the curvature 
as a measure of stability. The system becomes more stable if the absolute value of the 
curvature decreases. 

In the quantum mechanical regime (where rj >> 1) the curvature as a function of 
magnetic field has a maximum. This maximum occurs at stronger magnetic fields as the 
magnitude of fugacity is increased. Stability interpretation of curvature can not explain 
this maximum. 
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